The role of the Dirac sea in the quenching of Gamow-TeUer strength 
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The issue of Gamow-TeUer sum rule in atomic nuclei is reexamined in the framework of a rel- 
ativistic description of nuclei. It is shown that the total Gamow-TeUer strength in the nucleon 
sector is reduced by about 12% as compared to the prediction of a non-relativistic description. This 
reduction is due to the antinucleon degrees of freedom. 
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The quenching of the observed Gamow-TeUer (GT) 
strength in atomic nuclei is a long standing issue both 
experimentally and theoretically. On the experimental 
side the recent measurements 1] have pinned down the 
amount of strength reduction to about 10% in medium- 
heavy nuclei. Theoretically there are several possible 
sources of quenching ranging from multiconfiguration 
spreading^ to coupling to the A-hole sector|3]. The 
nuclear GT transitions are essential to many important 
processes in particle and astrophysics related to neutrino- 
nucleus interactions and therefore, further understanding 
of the GT quenching is needed. In this work we point 
out a new quenching mechanism due to the effects of the 
Dirac sea states. 

We start from the model-independent Ikeda-Fujii- 
Fujita sum rule0 



( IQ+Q-I >-( IQ-Q+I > = 2(iv-z), (1) 



where 



i 

with the isospin operators t± = (tx ± iTy) /V2. If we 
assume that there is no ground-state correlation, 



0+1 > -0, 

and Eq.J^I becomes simply 

( \Q+Q-\ )=2{N-Z). 



(3) 



(4) 



In this work we study the GT sum rule in a relativistic 
model where the nucleus is assumed to be composed of 
Dirac particles bound in Lorentz scalar and vector po- 
tentials. This model has been extensively studied for the 
last 30 years, and shown to describe very well both nu- 
clear reactions and structure phenomenologicallyj^. We 
will show that about 12% of the sum rule value Eq.QJ 
is taken by nucleon-antinucleon states with high excita- 
tion energy. As a result the strength of the giant GT 
resonance is quenched by this amount. Moreover, even 



if there is no ground-state correlations, Eq.||2Jl does not 
hold, since there are nucleon-antinucleon excitations. 
Let us start with the field: 
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where the two terms describe the particles and the an- 
tiparticles, respectively. The suffix a denotes the spin 
and isospin quantum numbers as Uq,(p) — Uct(p) I''') , and 
the positive and negative energy spinors are written as 



w<t(p) 
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with Ep = a/M*2 -\- p2 and ^ is the Pauli spinor. The 
Lorentz vector potential does not appear explicitly in the 
present discussions, while the Lorentz scalar potential 
does in the nucleon effective mass M*. Using the above 
field, an operator F can be expressed as 



F= / d^a;V'(x)/V'(x) 



(7) 



where / is considered to be some 4x4 matrix with no 
momentum-transfer dependence. From now on we will 
discuss nuclei where > fcp, and kp being the Fermi 
momentum of neutrons and protons, respectively. Then, 
we have for the t_ and t+ excitations 



F\ 



F+ 



V2 j (fp ^(^u„(p)rw^'(p) 4p(p)a^/„(p) 
+ u,{Y>)rv,,{^Y>)aUp)bln{-v)) I >, (8) 

V2 Jd^pY,^Ap)rv,'{~p) 

<7<7' 

xat„(p)foLp(-p)l ), (9) 
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where F stands for the 4x4 matrix other than the isospin 
operator. Straightforward calculations yield 



( \F-F+\ ) = - 



2V 
2V 



(27r)3 



(10) 



where we have used the fact that S{0) — V/{2Tr)^, V — 
A (Stt^ /2kp) being the volume of the system, and the 
abbreviations: 



with 



A± = Tr {FA+rA^ 



A+ = ^Ua{p)Ua{p) 



M* 



2Er, 



M* 



2E^ 



(11) 

(12) 
(13) 



and — {—Ep, p). The step function is also defined as 

^p"* = (^{ki — IpI) for i = p or n. In the case of the GT 
transition F — 757^, the traces of Ea. ((TTl) are calculated 

as 



A+ = 2 



A_ = - 2 
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(14) 



Finally, we obtain the matrix elements for the /3_ and /J^ 
GT transitions 



\F-F+ 



(27r)3 J E^ \ P 

X "0(")(Ar2+p2)+(p2_p2)] ^ (15) 

W fd^p 



(l-<))(p2-p^)(16) 



As seen from the above derivation, the second term in 
the square brackets of Ea. ((T3|l and the r.h.s. of Ea. l(TC|) 
come from the nucleon-antinucleon excitations, and each 
of them is divergent. However, the difference between 
Eq. ((T5)l and Eq. (fTB)l gives the sum rule value as it should 
be, 



\F+F-\) - {\F-F+\) = 2{N - Z). 



(17) 



If we neglect the antinucleons, then F+\ ) — while the 
first term in the square brackets in Ea. (|15|l gives 



(2^) 

A 
P" 



^ + 2Af*2 f fc„ - M* tan-i 
3 V M* 



(18) 



Expanding the above equation up to first order in (fcn 
kp), we obtain 



« 2 ( 1 - ) (A^- Z) , (19) 



where vp is the Fermi velocity 
fcp 



Ef 



Ev = ./M*^ + k^ 



and we have assumed as usual that 



(K - fcp)/fcF ~ 2{N - Z)/3A. 



(20) 



Thus, the sum rule value in the nucleon space is 
quenched. In most of the relativistic modelsQ the nu- 
cleon effective mass is about 0.6M, which yields vp = 
0.43 for kp = 1.36 fm~^. This corresponds to a quench- 
ing of 12%. If we use the free nucleon mass, the quenching 
is about 5%. The fraction of sum rule taken by nucleon- 
antinucleon excitations: 



[( \F+F^\ 



Jnn 



2^«l (N- 



Z) (21) 



must be distributed over the higher excitation-energy re- 
gion, but we need the renormalization for more details. 

The amount of quenching depends on the quasi- 
particle distribution, since the commutation relation tak- 
ing only the nucleon degrees of freedom becomes 



Pv 



El 



aL(p)a^n(p) - aTp(p)a„p(p)j(22) 

If quasi-particles have high- momentum components, the 
amount of quenching would be increased. In the limit 
fcp ^ 00, the ratio of the GT strength of nucleon excita- 
tions to that of antinucleon excitations is 1 to 2. 

We note that in the case of the sum rule for the isobaric 
analogue state, which is similar to Ea. HlTI) . there is no 
contribution from the Dirac sea. This fact is shown in 
the same way as for the GT sum rule. Taking P = 70 in 
Ea. (|ll(l . we obtain A+ = 2 and A_ = 0, which lead to 

{\F+F_\)^2{N~Z), {\F^F+\)^0 (23) 

instead of Eqs.lCHI) and 

Since a part of the GT strength is taken by the nucleon- 
antinucleon states, we investigate the effects of the Pauli 
blocking terms on the strength of the giant GT resonance 
in RPA. It has been shown that the Pauli blocking terms 
play an important role for some observables in the rela- 
tivistic model 0- We assume that the RPA correlations 
are induced through the Lagrangianfsj: 



ibTn , 9b = 9 



I J IT 
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where g' stands for the Landau-Migdal parameter. Al- 
though the way to introduce g' in the relativistic model is 
model-dependent we take this form which reproduces 
the non-relativistic result for the excitation energy of the 
giant GT resonance, as will be shown later. Then, the 
RPA correlation function TTrpa for the r_ excitations is 
written in terms of the mean field one 77, 

^rRPA(r^ , r^) - [c/- Y^^(A+ , r^), (24) 

where U denotes the dimesic function: 

u^"" = 9'"'-X5nir'^,n) (25) 
I 



with X5 = 55/(27r)^. The general form of the mean field 
correlation function including the Pauli blocking terms 
can be found in Refs.0,0. We follow the usual proce- 
dure of taking into account all the terms containing the 
density-dependent part Gd of the single-particle propa- 
gator G, but neglecting terms like rnGyFi^jGy (Gy is the 
Feynman part of G) which may be divergent |9|, [lOj and 
should be treated by a renormalization procedure. For 
n{r^,r'^ ) one obtains 



^(PO - Ep) 



{p + q? 



M 



*2 
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(p-g)2-M*2 + i£ P 



Epl^p+q 



(26) 



where t^'^ {p, q) is 

t^''ip,q) = A(^g^''{M*'+p'+p-q) 



(27) 



The above equations show that the dimesic function can 
be separated into transverse and longitudinal parts. In 
the case of q = 0, the determinant can be written as 

det U = - iDrf Di^, (28) 

where we have defined 

Dt = i + x^n{rl,rl), (29) 
Di. = i-x^n{rl,rl). (30) 

The eigenvalues of the excitation energies are obtained 
from the real part of the dimesic function, 

det Re C/ = 0. 

The longitudinal part has no solution, while the trans- 
verse part yields for go < M* with Eq. 1(20)) ^ 



WGT 



l-2w|/3 ^klN-Z 



1 -55P(fcF)'^^37r2 2 A ' 
where P{k-p) stems from the Pauli blocking terms. 



(31) 



P(M = 3^i?l(^-F~«|-^(l-«i)logi + "^ 



Tr^K4 1 + + • 



(32) 



In the non-relativistic limit w| 
gives 



^ 1, the above equation 



37r2 2A 



(33) 



which is just the excitation energy of the giant GT res- 
onance obtained in the non-relativistic modeljllj. If we 
use the value g'— 0.6 12], the Pauli blocking effect on the 
excitation energy is less than 0.5%. 

The strength of the giant GT resonance is given by the 
imaginary part of the RPA correlation function. Writing 
n{qo) — n{r^,r^ ), we expand the transverse dimesic 
function £'t('7o) at qg = wgt, 



Driqa) = ax (go - ^gt) + 



dDT 



dqo 



Then, the RPA correlation function is 

1 niqo) 



^rpa('7o) 



ax qo — ^GT 
which gives the imaginary part 

TT 



Imi7RPA((7o) 



S{qo ~ loot)- 



(34) 



(35) 



(36) 



The relationship of the imaginary part to the 
strength 10], together with Ea. (|31|l . provides us with the 
strength of the giant GT resonance: 



GT 



{l^g5P{kF)r 



1--v'f]{N~Z). (37) 



The above result shows that S'gt is almost equal to the 
one in the nucleon space given in Ea. H19|) . since P{kp) 
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is negligible. Thus, the Pauli blocking terms yield al- 
most no effect on the strength of the giant GT resonance. 
When we neglect P{kp), S'gt is independent of 55 and it 
exhausts the sum value in the nucleon space. 

It should be noted that the Pauli blocking terms 
change the sum value of the unperturbed states. The 
imaginary part of the mean field correlation function is 
given by F,q.^^ at q = 0: 



imn{rl,r^ ) 

— TT / d p ■ 



+ t^\p,-q)S{q^o-2pm)e^J'^ 



(38) 



When we employ the relation: 



S{ql ± 2^090) 



1 

2?^ 



{diqo)+S{qo±2po)) , (39) 



the first and second ^-functions of the r.h.s. correspond 
to the sum values from the nucleon space and from the 
Pauli blocking NN-excitations, respectively. The latter 
yields the sum value for the /3_ transition: 



'5'nn 
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0{V) 



(40) 



which is negative, and equal to the density-dependent 
part from NN-excitations in Ea. ll5() . 

In conclusion, we have shown that, in a rclativistic 
description the total GT strength in the nucleon sector 
is quenched by about 12% owing to the antinuclcon de- 
grees of freedom. The effects on the giant GT resonance 
energy are negligible. Of course, the strength in the nu- 
cleon sector can be further spread out to the tail of the 
giant resonance by configuration mixing effects similarly 
to the non-relativistic picture. It has been discussed for 
a long time that the GT strength may be quenched be- 
cause of the coupling of the particle-hole states with A- 
hole states 12]. The present discussions predict another 
source of the quenching. The recent experiments have ob- 
served about 90% of the classical sum rule value[3|. One 
should study in more detail how A-hole and NN states 
contribute to the observed quenching. It would be inter- 
esting to extend the present investigation to include the 



diverging FaG-prpG-p terms by a renormalization proce- 
dure for a consistent treatment in the future. 



Finally, we note that the Coulomb sum values are 
also strongly quenched by the antinucleon degrees of 
freedom 13]. This fact is consistent with recent analysis 
of experimental data 14]. Sum rules require the complete 
set of nuclear wave functions. This implies generally that 
antinucleon degrees of freedom have a possibility to mod- 
ify the sum rule values as compared to non-relativistic 
models. 
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